Magnetization of a D.C. biased quantum dot 
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Using a quantum generalization of the Onsager principle of microscopic reversibility, the magneti- 
zation of a system in a non-equilibrium steady state quantum dot is formulated as a response of 
the interaction energy to an external field. This formulation permits a direct and compact compu- 
tation of the steady-state magnetization of a non-equilibrium quantum dot as a differential of the 
interaction energy. Unlike the direct computation of the magnetization using perturbative Keldysh 
methods, this approach does not require the use of a point splitting procedure. Our results never- 
theless support earlier calculations made in the limit of zero field, and they support the survival of 
strong coupling to arbitrarily large voltages, both at zero field, and under the conditions where the 
chemical potential difference eV becomes equal to the spin-fiip energy in a field eV = q^lbB. 

PACS numbers: 



I. INTRODUCTION 

Quantum dotsS offer a unique opportunity to study 
the behavior of a simple, yet non-trivial strongly cor- 
related electron system out of equilibriumBllj In the 
Coulomb blockade regime, quantum dots with odd num- 
bers of electrons develop a magnetic moment that in- 
teracts non-perturbatively with the leads at low temper- 
atures. The unbiased quantum dot provides a classic 
example of strongly correlated electron physics: at low 
temperatures, the residual interaction of the dot spin with 
electrons in the leads gives rise to a Kondo effect, whereby 
the formation of a many -body resonance at the Fermi 
energy drives the conductance in the Coulomb blockade 
regime back up to the unitary limit. The response of the 
Kondo effect to a D.C. bias is a matter of some interest, 
and provides us with a chance to examine how a non- 
trivial strongly correlated electron system responds in a 
non-equilibrium setting. 

One of the key properties of the Kondo effect in equilib-, 
rium is the presence of a "running coupling constant"til 
whereby the antiferromagnetic coupling between the dot 
spin and the leads g{h) grows progressively as the en- 
ergy scale A is reduced. The Kondo effect exhibits the 
phenomenon of confinement. At high energies, the local 
moment is asymptotically free, weakly interacting with 
the surrounding environment, but at low energies the 
spin of the dot is "confined" : screened by the lead elec- 
tron. A single scale, the "Kondo temperature" Tk^ gov- 
erns the low temperature properties; Tk = D^/g e~^/^^ , 
where g is the 'bare coupling' between the spin of the 
dot and the leads and D ^ Tk is the electron band- 
width. Thus, for example, the magnetization at tem- 
perature T and magnetic field B is a universal func- 
tion M = m{T/TK, B/Tk), where M has a perturbative 
"weak coupling" expansion in g only when T,B ^ Tk ■ 
(See Fig. 1(a).) 



What happens to the Kondo effect when a current 
flows through a D.C. biased quantum dot? Some have 
suggested that the passage of electrons through the 
quantvHH |d«t will "decohere" the physics of the Kondo 
effectEHl3il3, dephasing the coherent spin-flip processes 
necessary for the Kondo screening process. The rate at 
which electrons pass through the quantum dot is given 
by I/e, where / is the current through the dot. The 
formation of a Kondo resonance requires that quantum 
processes are coherent on a time scale ~ Ti/Tk- Sup- 
pose that quantum processes on time-scales in excess of 
Tc ^ ejl are dephased by the passage of electrons, then 
when Tc ~ Ti/Tk, or when 



sTk 



the coherent formation of a Kondo resonance is expected 
to break-down. This would imply that the Kondo effect 
will break down at voltages in excess of the Kondo tem- 
perature. In such a picture, the magnetization of the 
quantum dot 

M = M{T/Tk, H/Tk, cV/Tk) 

should become perturbative in the coupling to the leads 
at arbitrarily low temperatures once eV » Tk- 

An alternative picture of the quantum dot argues that 
the Kondo effect is modifled, but not|-deahased by the 
passage of electrons through the dotHEj. According 
to this picture, although interlead Kondo processes are 
cut-off by a finite voltage, intralead Kondo renormaliza- 
tion effects continue unabated in each lead until the sys- 
tem enters a non-equilibrium strong-coupling regime at a 
renormalized temperature T^{V). If true, the DC biased 
quantum dot should display a Curie like susceptibility 
characteristic of a decoupled local moment even in the 
presence of a curent, until the temperature approaches 
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the renormalized Kondo temperature T^(V). 
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FIG. 1. (a) Field dependence of the crossover temperature 
T^{B) separating weak- and strong-coupling regimes of the 
equilibrium Kondo model. At absolute zero, for fields larger 
than the Kondo temperature the quantum dot re-enters weak 
coupling, (b) Voltage dependence of the crossover temper- 
ature T^^y), from earlier results of Coleman, Hooley and 
ParcoUetta. Note that goes to zero only in the limit 
V — > oo, i.e. for low enough temperatures, the model reaches 
a strong coupling state irrespective of the magnitude of the 
voltage. 

Such a possibihty rests on the observation that the Kondo 
screening process partially delocalizes the spin of the 
quantum dot into the leads to form a kind of "meso-spin" 

S*iA) = [/(A)5C/t(A), 

whose composite structure depends on the scale on which 
it is observed. Here, the idea is that same unitary trans- 
formation /7(A) which integrates out the high energy 
electrons also admixes and hence delocalizes the bare dot 
spin S into the leads. (Fig. 2.) In this picture, the rapid 
fluctuations of the bare quantum dot spin 5 at a rate 
comparable with r~^, are merely internal redistributions 
of spin within the scale dependent meso-spin, 

5| ^ {high energy singlet} -I- {e, h}p 



(s*(A)) = 



+ 



FIG. 2. Renormalized "meso-spin" which forms around 
a magnetic ion as high-energy electrons are systematically 
integrated out of the spectrum. The original spin of the bare 
local moment is partially delocalized into the conduction sea. 
Fluctuations of the bare moment at a rate lead to an 
internal redistribution of spin within the meso-spin, but do 
not cause the composite object to flip. 



If this picture is correct, the magnetization of the DC bi- 
ased quantum dot is expected to remain perturbatively 
close to a Curie law at large voltages down to a renor- 
malized Kondo temperature T^{V), where a new kind of 
non-equilibrium Kondo effect will take place. 

These discussions suggest that a careful computation 
of the quantum dot magnetization may provide a key to 
understanding whether a Kondo effect takes place in the 
quantum dot at a large voltage bias. In this paper, we 
present an extended discussion of earlier perturbative cal- 
culations of the quantttm dot magnetization by Coleman, 
Hooley and Parcollett3 that suggest that strong coupling 
physics extends to arbitrarily high voltages in the D.C. 
biased Kondo model (Fig. 1 (b)). Two of the thr^e orig- 
inal authors of that work, ParcoUet and HooleyE3 have 
recently argued that the key results are incorrect, claim- 
ing that the perturbative expansion of the magnetization 
breaks down at zeroth order in the coupling constants 
when a current flows through the quantum dot. Part of 
the reason for doubts about the earlier work stemmed 
from the necessity of introducing a "point-splitting" pro- 
cedure to control divergences in the perturbation theory 
associated with the conservation of the total magneti- 
zation. In this paper we present new arguments in fa- 
vor of the-|Original conclusions of Coleman, Hooley and 
ParcoUetllj. Our arguments avoid some of the earlier 
difhculties associated with point-splitting by taking ad- 
vantage of a quantum mechanical extejusion of Onsager's 
principle of microscopic reversibilityndH to steady-state 
non-equilibrium quantum mechanics. 



II. MODEL AND KEY RESULTS 

Under conditions where the quantum dot has almost 
integral valence, interactions between the local moment 
on the dot and its leads can be described by a Kondo 
model, of the following form: 

H £kcl^l^^C,nka +Hl - B MtOT, 



ikCT 



Hi ^ J'0'^Q(O)(Ta/3V'/3 • S. 



(1) 



Here, cj^j.^. creates an electron in lead m e {L, R} with 
momentum k and spin cr, and J is a positive (antiferro- 
magnetic) Kondo coupling constant derived from virtual 
charge fluctuations on and off the quantum dot. For sim- 
plicity the spectrum of the leads on either side of the dot 
is taken to be identical and 



V'a(O) 



fc 



(2) 



describes a linear combination of electrons in the left and 
right hand- leads with aj^ + aj^ = 1. In equilibrium the 
left and right-hand leads are filled to the same chemi- 
cal potential /i, but out-of-equilibrium, the left and right 
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hand leads differ by an amount equal to the driving volt- 
age 



(3) 



Finally, note that we have included the magnetic field 
term in the interaction part of the Hamiltonian, where 

MtOT ~ 2Sz + Mieads 

is the total magnetization of the leads plus dot, 
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is the magnetization of the leads and n^kcr = c^Ako-CAko- 
are the number operators for electrons in the leads. 

The interaction term in this Hamiltonian can divided 
up into inter and intra-lead terms, as follows 



(4) 



where 



E 



ma^crcr'Cm'a'j ■ S + (R ^ L) 



k,k' ,a".cr' 
H trans — JlR. > I C 



/ . \ ^RWa'- 
k,k',cr,a"' 



^aa'CD^'a' ) ' S + {R ^ L) 



describe the spin-flip reflection and transmission of elec- 
trons between the leads and the coupling constants 



Jl 
Jrl 



Jr. = oij^J 



Jrl — OLRahJ 



(5) 



It is sometimes useful to consider the coupling constants 
Jl, Jr and Jb.l as independent variables. 

In equilibrium, the Kondo model has been solved ex- 
haustively by a number of different methods. Central 
to our whole understanding of the Kondo model, is the 
renormalization group. Unfortunately, a complete frame- 
work for the renormalization group out-of-equilibrium is 
not yet available. Much of the current understanding 
of the Kondo model out-of-equilibrium is based on re- 
summed perturbation theory. A second approach to the 
problem has been to use a strong-coupling analysis. Each 
of these approaches is faced with its own particular bias. 
In a resummed perturbation theory, certain classes of 
processes are selected, and there is genuine danger that 
cancellations inherent to the full problem are lost in the 
process of selecting diagrams. By contrast, the strong 
coupling approach relies heavily on the linearization of 
the band, and the assumption that the pertinent band- 
electron cut-offs in strong coupling are larger than the 
applied voltage. However, we know from the equilibrium 
Kondo model, that in some sense, the band-width of the 
strong-coupling problem is given by the Kondo temper- 
ature, so it is not clear how this last assumption can be 



reconciled with a voltage larger than the Kondo temper- 
ature. 

In the face of these uncertainties, we argue that the 
only unbiased way to address the issue of whether the 
DC biased Kondo model enters a strong-coupling regime, 
is to revisit perturbation theory in a non-equilibrium 
context. There are two obvious physical variables that 
can be used to characterize the Kondo effect out-of- 
equilibrium: the current through the dot, and the tp=. 
tal magnetization. Perturbative studies of the currentE2l 
through a quantum dot indicate that the expansion 

I^'-f{T/TK,eV/TK) 

h 

enters a weak coupling, perturbative regime when eV » 
Tk- At first sight, this might be taken as evidence that 
the Kondo effect returns to weak coupling at large voltage 
bias. However, the current between two leads at large 
voltage bias eV » Tk involves electrons far from the 
Fermi surface of each lead, and these electrons would be 
weakly coupled to the dot spin even under equilibrium 
conditions. 

We shall argue that a more useful quantity to charac- 
terize the behavior of a quantum dot out of equilibrium, 
is the static magnetization. We shall take the magnetiza- 
tion M to refer to the impurity contribution to the total 
magnetization, defined as the difference 

M = {MtOt)j - {Mieads) J =Q- 

In the region of weak coupling max{T, H)/Tk » 1, 
the equilibrium magnetization has a perturbative expan- 
sion. For low fields, M{T) = x{T)B, where 



xiT) 



1 

T 



l-2Jp-A{Jpf 



r)e(3/4+7) 

\og(— ) 

27rr ' 



(6) 



where 7 = 0.5772 = is Euler's constant, D is the 

band-width and p is the density of states in each lead. 
In the same cut-off scheme, the high-field magnetization 
B » T \s given by 



M = 1 



l-Jp-2{JpfH^)- 



(7) 



The cross-over to strong-coupling is defined by the scale 
at which the second and first order corrections to the 
magnetic susceptibility and magnetization become equal. 
This defines a field scale 



2B = Tk = Deyij>[~l/{2Jp)] 
and a temperature scale. 



To^D- 



27r 



■exp[-l/(2Jp)] ^ QMTk 



(8) 



(9) 



In this paper, we generalize these results to the case 
of finite DC bias, computing the magnetization to sec- 
ond order in the coupling constant at finite temperature. 
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field and voltage. In order to do this, we need to define 
what we mean by the magnetization. We consider the 
effect of turning on the the Kondo coupling constant and 
computing how the magnetization changes in the steady 
state, as we slowly increase the coupling constant up to 
its maximum value. We shall argue that under rather 
general conditions, the change in the magnetization due 
to the Kondo interaction can be written in the form 



AM 



_d_ 

"dB 



— {'Hl)j',H 



where {Tii) j^h is the steady state value of the interaction 
energy, evaluated at Kondo coupling constant J, in a 
magnetic field B. The advantage of this formulation, 
is that it permits a ready generalization of equilibrium 
methods to a non-equilibrium, steady state situation. 




FIG. 3. Showing the region of strong coupling defined by 
the equality of the first and second order terms in the ex- 
pansion of the magnetization. Strong coupling persists to 
arbitrarily high voltage at B = and for gfisB ~ eV. 



The main result of our calculation is succinctly cap- 
tured by the zero temperature magnetization, which for 
the symmetric case Jr = Jl = Jrl = J/^, takes the 
form 



M 



l-Jp- {Jpf In 



2B^\{eV)^-{2By\ 
This expression is perturbative so long as 



2B^\ieVr~{2Br\>{Ti 



K 



(11) 



which defines the region shown in Fig. 2. 

As can be seen from Fig. 2., large voltage identifies 
two regions of strong-coupling: 

• Zero field Kondo effect 2B < {TxY/eV, corre- 
sponding to a strong coupling region associated 
with the intra-lead Kondo processes that are not 
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cut-off by a finite voltage. Earlier workia^t^ sug- 
gested the corresponding fixed point is a two- 
channel Kondo model where 

More recent work by Rosch et aS and Schiller et 
alEZl has suggested that the entry into the strong- 
coupling two-channel regime is ultimately cut-off by 
decoherence effects. We shall return to this point 
in the discussion. 

• Finite field Kondo effect. This occurs around 
- 2B = gpsB, \eV - 2B\ < Tf^/{eVf, and 
also extends up to arbitrarily large voltage. Here 
the Kondo effect involves the degeneracy between 
a spin flip of the quantum dot and the transfer of 
an electron between the fermi surface of the two 
leads. The effective Kondo model for this situation 
presumeably resembles a one-channel Kondo model 
at zero voltage bias, where the down electrons on 
high voltage lead and the "up" electrons on the low 
voltage lead can scatter elastically between leads 
by exchanging spin with the local moment. We ex- 
pect that the physics at energy scales lower than 
A max(f; /i^i?, eV), will involve only these elec- 
trons, and the effective model will be a one channel 
Kondo model of the form 



(12) 



where the spin S sits in a zero field, and the chemi- 
cal potential difference between the L | and R t 
electrons is A/i* — A/i — gpsB. In this situa- 
tion, we expect the spin on the quantum dot to 
be perfectly quenched at low temperatures . Small 
changes in the voltage will affect magnetization, 
however, without a mechanism for the transfer of 
electrons from one lead to another, this finite field 
Kondo effect will not, it appears, lead to a unitary 
differential conductance. 

At finite temperatures, our general result takes the 
form 



(10) M ={l~ J) tanh {B[l~ J\/T) 



2[{JRf + {JLf]^ 



dB 
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where ~ J^p, J = {Jr. + Jl) and 

1 
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Re / dx 



(2cosh(f))2' 



— 27ri(T Inf (x + (T& + 7w) 



(T,7 — ±1 
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order in perturbation theory, is a signature that the intra- 
lead Kondo effect continues unabated at temperatures or 
fields smaller than the voltage. 

The corresponding boundary where weak-coupling per- 
turbation theory breaks down is shown in Fig. 4. 
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FIG. 4. Illustrating the Kondo effect at a finite field 
qHbB = eVo- At absolute zero, we expect the ground-state to 
be a singlet and small departures in the voltage from Vb will 
then lead to a magnetic polarization of the local moment. 

In the zero field, finite temperature limit, this general 
result reverts to the form M — xB, where 



1 

T 



-S\Jlb.\ 




where 7 = 0.5772 . . . ~ — V'(l) is the Eulcr constant. The 
crossover function (l){x), in terms of digamma functions 

ip{z), is 



= Re 



dy 



-co 4cosh2(f) 



^iy + x)-ij{y) , (15) 



where i'{x) = ^p{^ + i-§^)- 

The second-order terms in describe the leading 
logarithmic enhancement of the Kondo coupling. Terms 
of order J^^j involve inter-lead processes and, as ex- 
pected, the logarithmic divergence in these terms is cut 
by the voltage. (To see this, note that —(l){V/T) ^ 
— hi{V/T) for y ^ T which cancels the logarithmic tem- 
perature divergence.) By contrast, the intra-lead terms of 
order and J| are completely unaffected by the voltage 
V, which guarantees that the leading logarithmic diver- 
gence survives at arbitrarily high voltage. This is easily 
seen in the large-y form of the susceptibility. 



X. = - ^ 1 - 2( J« + Jl) - 4 {Jl + Jj, 



-8|Jiflrin( I 



(17) 



The survival of the leading logarithms in the suscepti- 
bility at arbitrarily high voltage is consistent with the 
zero-temperature, finite field results, and at least at this 
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Fig. 4. Illustrating region of strong coupling for finite 
field, temperature and voltage. 



III. MAGNETIZATION AS A RESPONSE 
FUNCTION 

Since we wish to consider the case where the voltage 
Vsd takes an arbitrarily large value, linear response meth- 
ods are not appropriate; instead, we must resort to the 
Schwinger-Keldysh formalism for non-equilibrium field 
theory. The SK mcthodcj divides the Hamiltonian of 
the system into two parts: 



H = Ho + Hi + X{t)A, 



(18) 



where, for convenience we have introduced a source term 
with field X{t) coupled to quantity A. It is assumed that 
we know the exact Green's functions for the operators 
under Hq- The interaction Hi and source term are re- 
garded as perturbations. 

In the SK approach, the expectation value of an op- 
erator A{to) in the interaction representation A{to) 
eiHot^Ae- " 



IS 



{A{to)) = {P A{to)e 



J K 



dt 



H,+\{t)A{t) 



(19) 



where K denotes the two branch Keldysh contour run- 
ning from —00 to -l-oo back to t = —00 , (Fig. ^) P de- 
notes path ordering along this path and (...)= Tr[po • ■ •] 
is average with respect to the density matrix po ( in the 
absence of the interaction) in the distant past. 
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FIG. 5. The Keldysh contour. The closed contour ehm- 
inates the need to normahze the path-ordered exponential. 
Fields are defined on both the upper {tp{t+)) and lower 
{tp{t-)) contour. Along the path, the time t- is taken to 
be "later" than the time t+. 

In the conventional SK approach, the expectation 
value of some quantity is computed directly, using di- 
agrammatic methods. We may regard this as a direct 
"measurement" of the quantity in question. In this paper 
we simplify our computational methods by implement- 
ing a new approach which regards expectation values as 
a response to external fields. The inspiration to go from 
a "measurement" to a "response" based approach is de- 
rived from equilibrium quantum mechanics. 

A major difficulty in the Keldysh approach, is that the 
measurement vertex is much more awkward to use than 
the "response" vertex. For example, in the Larkin Ov- 
chinikov basis the response vertex is diagonal whereas the 
measurement vertex is off-diagonal (— iri). This leads to 
serious difficulties when dealing with a conserved order 
parameter, for double poles which appear in the pertur- 
bation theory can not be absorbed using integration by 
parts. The same double poles are easily handled by re- 
garding the expectation value of a quantity as the re- 
sponse to its conjugate field, because the response vertex 
is proportional to the unit Keldysh matrix, and between 
two bare Green's functions it satisfies the relation 

gg = -d^g. 

This permits us to eliminate double poles in the per- 
turbation theory using integration by parts. It is also 
extremely useful in the development of Ward identities. 

In equilibrium thermodynamics, we first calculate the 
Free energy, and then we take the derivative to compute 
the magnetization 



M ^ - 



dF 
dB 



Can an analogous approach be used out-of-equilibrium? 

We now make a digression on the nature of the non- 
equilibrium steady state. One of the key features of ther- 
mal equilibrium, is the notion that the equilibrium steady 
state does not depend on its past history. In particular, 
if we turn on some interaction of strength g, Hj = ghj, 
where hj is a dimensionless interaction operator, and we 
couple internal degrees of freedom Ai (z = f,n) to cor- 
responding external fields Ai, then the amount of work 
done on the system 



AW 



J {hi{t))dg{t) + (A,(t))dA„ 



where we have used the summation convention on i, does 
not depend on the detailed path P over which the cou- 
pling constants {g, A;) are adiabatically incremented to 
their final value. This implies that the "curl" defined by 
the following functional derivatives must vanish 



5{hi{t)) 5{Mt)) 

5K{t') 5g{t') 

5{A,(t)) 5{Mt)) 

5\,{t') 5\j 







(20) 



^From linear response theory, we can relate these func- 
tional derivatives to the corresponding response func- 
tions, 

^J^^^-^{[hJ{t),A,{t')])6{t-t') 
^^^^-^{[Adi)M{t')mt-t') 
^J^A=-^(^[A,{t),A{t')])eit-t') (21) 



from which it follows that 



^i{[hi{i),M{2)])e{i 

-z([A,(I),A,(2)])e(f 



-i{[A,{l),hi{2)])e{l~2), 
-i{[A,{l),A,{2)])e{l-2). (22) 



These relations are the quantum mechanical iCopiiterpart 
of the famous Onsager reciprocity relationsEiltJ. They 
reflect the microscopic reversibility of the equations of 
motion and the absence of any "arrow of time" in thermal 
equilibrium. 

We now argue that Onsager's reciprocity relation will 
continue to hold for an important sub-class of variables 
in the non-equilibrium steady state. For our discussion, 
we consider a quantum dot that is coupled to two very 
large baths of electrons ("leads") at different chemical 
potentials /il aud /i^ where ijll > Mfi The entire system 
is isolated. Suppose we connect the quantum dot to the 
leads at time i = 0, then we expect that after some equi- 
libration time Ti , which is basically independent of the 
size of the leads, the system will arrive at a steady state 
where a current of electrons flows from the left, to the 
right-hand lead. This steady state will persist for a long 
time T2{L) until a substantial fraction of the additional 
electrons on the left lead have flowed into the right lead. 
Typically, the time T2{L) over which the steady state will 
persist, will diverge as L ^ oo. The steady state value 
of some variable A is then given by 

{A) = hm {A{t)) 

L — >QO 

with the understanding that T2(L) >> t » ti. 

The steady state involves a persistent current from left 
to right, and evidently involves an "arrow of time". Nev- 
ertheless, we expect that the nature of the steady state 
does not depend on the history of how it is arrived at. In 
particular, suppose we the steady state is arrived by adi- 
abatically turning on an interaction Hj = ghj between 
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the leads, and by coupling source terms \j to various 
quantities Aj associated with the dot. Since the system 
is closed, when we adiabatically change these variables 
we can define the amount of work done in reaching the 
steady state 

^Wne = I {hi(t))dg(t) + {Adt))dX,, 



We emphasize that since the system is completely iso- 
lated, the amount of work done in reaching the final state 
is literally the change in the total energy of the system: 
a completely well-defined quantity. 
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FIG. 6. Illustrating a closed quantum mechanical system 
in which the non-equihbrium steady state of a quantum dot 
can be studied, (a) Quantum dot coupled to two isolated 
leads, initial in thermal equilibrium at two different chemical 
potentials, via an interaction term Hi — g{t)hi. The coupling 
constant g (t)is switched on at time t = 0. (b) After an initial 
transient of time ti, the current reaches a steady state value 
which only depends on the chemical potential difference and 
the other fields (e.g the magnetic field) applied to the dot. 
This steady state persists for a large time T2{L) which goes 
to infinity as the size of the leads is sent to infinity. 

If the steady state is independent of the path by which 
it is arrive at, then it is not unreasonable to expect that 
there is a sub-class of variables {Ai} for which the work 
AWne is path independent. If this case, that clearly the 
Onsager relations given above must extend to the non- 
equilibrium steady state for this subclass of variables. Of 
course we do not expect the reciprocity relation to extend 
to all variables, as it docs in thermal equilibrium, because 
this would mean that the arrow of time is completely in- 
visible. However, the idea that there exists a subclass 
of "conforming variables" whose correlation and reponse 
functions are completely symmetric in time, if true, could 
be invaluable for studying the non-equilibrium steady 
state. 

This reasoning has motivated us to, recently propose a 
"Quantum Reciprocity Conjecture"E3 : 



In the non equilibrium steady state, the set 
of quantum mechanical observables contains 
a non-trivial subset V of "protected" quan- 
tum observables V — {ai, 02 ... , a„} whose 
correlation functions in the steady state are 
insensitive to the arrow of time, and which 
consequently satisfy a quantum mechanical 
analog of the Onsager reciprocity relations 

(Kl),5(2)]) = ([6(l),a(2)]), ia,b eV). 

Consider the retarded and advanced Green functions be- 
tween protected variables. 



Gif^^=T^([«(l),6(2)])^±(t2-ti) 



(23) 



where 0±{t) — 9{±t). Since a and b are hermitian, these 
are real functions {G^"^{t) = [G^'^{t)]*). The conjec- 
tured Onsager relations mean that in the steady state, 
they also satisfy 



'ah 



(24) 



where the order of the subscripts and time variables is 
important. If we write = G^* in the first relation, 
and then Fourier transform, we obtain the more familiar 
result 

GU^) = G^,(a;)* 

which means that the retarded and advanced Green func- 
tions of protected variables share the same spectral de- 
composition 



Gir^(-) 



dE 



1 



TT uj — E ± i6 



AabiE) 



where Aab{E) = ±Im[Gi^'^) (£;)]. 

Onsager used his reciprocity relations to propose 
that an effective "Free .eaesgy" might exist in a non- 
equilibrium steady stateEllEHl, an idea_which has its ori- 
gin in early works by Lord Rayleigho. Recently, Der- 
rida, Lebowitz and Speer have shown that this idea holds 
rigorously in certain simplified non-equilibrium classical 
modelsE3. Provided that the set of protected quantum 
variables includes the interaction hj, then it must fol- 
low that one can define an effective Free energy from the 
virtual work done in reaching the steady state 



AF[g,{Xj}] = {hi)dg + {A,)dX, 
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^dg + {A,)dX, 



(25) 



We have tested this conjecture, on a simple, non- 
interacting resonant level modelEj, where the hybdiza- 
tion with the leads Hj = ghj, the occupancy, magne- 
tization and spin current through the dot are found to 
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satisfy mutual Onsagcr reciprocity relations in the pres- 
ence of a finite steady state charge current. The ciirrent 
operator does not lie within the class of protected vari- 
ables, nor indeed does any operator which changes the 
balance of coupling between the leads, or with a separate 
heat bath. Nevertheless, the existence of a finite sub- 
class of protected variables has been confirmed in this 
non-interacting example. 

We now apply the above arguments on the assumption 
that the interactions inside the quantum dot are also pro- 
tected variables. If this is the case, then interactions can 
be adiabatically turned on inside the quantum dot, whilst 
preserving the reciprocity relations. This assumption en- 
ables us to extend the classic equilibrium expression for 
the change in magnetization as a consequence of coupling 
to the leads. 



AM : 



d_ 
dB 



I 



dg 



{Hi) 



(Egn) 



to the non-equilibrium situation 



_d_ 
dB 



I 



This "energetically defined" magnetization quantity mea- 
sures the response of the interaction energy to the exter- 
nal field in a the non-equilibrium steady state and it is an 
interesting quantity in its own right. We shall now show 
that this quantity is directly related to the change in the 
non-equilibrium magnetization due to interactions, 



AM 



NE 



AM. 



provided Quantum Reciprocity holds between the inter- 
action and the magnetization. 

Consider a non-equilibrium system in the steady state. 
To be careful, let us focus on some variable in the Hamil- 
tonian, call it A, which is conjugate to variable X{t), so 
that the Hamiltonian has the form H = Ho + Hj + 
\{t)A{t), where we implicitly assume that the interac- 
tion Hi = ghj has strength g, which we ultimately set to 
unity. Imagine turning on the interaction, the change in 
the A is then given by 



dAjto) 
dg 



{[A{to),Hj{t')])dt' 



where the response function is to be evaluated at non- 
zero g. If quantum reciprocity holds, 

-i{[A{to),Hj{t')]) = -i{[Hi{to),A{t')]), 

so that we can now put 

= -^- f" {[Hjito), Ait')])dt' (26) 
^9 9 J —oo 

But we can simplify this expression by writing 



-i{[Hi{to),A{t')])e{to-t') = 



S{Hi{to)) 
SX{t') 



so that 



dA{to) 



9 J-, 



dt 



MHi{to)) 



(27) 



dg g J-oo 
^ l d{Hi{to)) 

g d\ ■ 
Finally, integrating over g^ we have 

We will apply this expression to the magnetization in 
the non-equilibrium Kondo model, making the replace- 
ment A^ M and A — > —B and g — > J to obtain 

This expression deserves some discussion. In deriving it, 
we have assumed that the magnetization can relax to its 
new steady state value : a very delicate assumption! 

The magnetization is actually a conserved quantity 
which commutes with the full Hamiltonian [M, if] = 0, 
and so strictly speaking, by coupling the dot to the leads, 
we fail to relax the total magnetization. To avoid this 
paradox, we must surreptitiously introduce a regulated 
magnetization operator whose matrix elements lie be- 
tween states which are not degenerate, as in the case of 
a conserved operator, but whose energies are split by an 
energy scale s. Such an operator will consequently relax 
on a time scale = %/ s. It is important that this time- 
scale is faster than the adiabatic time scale ^ fi/S 
at which the interaction is switched on, where 5 is the 
regulator used inside the Green functions, i.e. 

Tr = ?i/s << %/5 = Ta- 

This is done in practice by weakly connecting the leads 
and the local moment to a heat bath, e.g. by consider- 
ing the local moment to weakly coupled to a spin chain. 
Formally, at the end of the calculation, we take 5 ^ at 
fixed s, then let s go to zero: 

M = lim lim(M)s s 

Formally, a precisely similar regulation is required in 
equilibrium. For example, when we consider the Pauli 
spin susceptibility, which is computed by taking the spin 
susceptibility at a finite q which is ultimately set to zero. 

XPauH = lim lim x(9,w) 

Fortunately, we need only consider this process in the ab- 
stract, because the computation of the interaction energy 
in a field does not present one with any divergent singu- 
larities. DifHculties in earlier calculations associated with 
point-splitting are completely eliminated. 
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IV. DEFINITION OF GREEN FUNCTIONS 

Let us denote by t+ a time lying on the upper branch of 
the contour, and by t_ a time lying on the lower branch. 
Then we may formally define four Green's functions: 



G++it,t') 
G+-{t,t') 

G—{t,t') 



-z(r^(<+)Vt(t;)); 



Here, T is the usual time ordering operator, while T is 
the anti-time-ordcring operator; we consider the opera- 
tors V' to be fermionic. The Green's functions in ( p8| ) are 
not all independent. In fact, there arc only two indepen- 
dent Green's functions here. This is most easily seen by 
employing a transformation into the so-called "Larkin- 
Ovchinnikov" basisEj. 

Let us define the matrix G', whose entries are the four 
Green's functions of (OS 



G' 



G+- 
G— 



(28) 



The transformation to the Larkin-Ovchinnikov basis is 
defined as 



G' ^ G = UG'V, 

where the matrices U and V are given by 



U = 



1 I 



V2 VI 



V = 



1 



x/2 



1 1 

-1 1 



(29) 



(30) 



The elements of this matrix Green- function correspond to 
the Green-functions of the symmetric and antisymmetric 
fields 



Defining 
then 

G(l,2) = 



-^(^(t+)-^(t_)). 



(31) 



^-(0 =^^(t)Ti = (V-^A^^s) (32) 



- c,/,t 



-i(P*(l)^'(2)) 



-i(P 



(33) 



where P denotes "path-ordering" of the fields along the 
Keldysh contour. It is easily shown that the joint fluc- 
tuations in the antisymmetric fields identically vanish, 
(Appendix A) 



and that the remaining matrix elements can be written 
as follows 



^ " I G^ 



(35) 



where the physical Green's functions are given by 
G''it,t')^~t0{t-t'){{m,i^Ht')}); 

G^{t,t')^29{t'^t){{m,i^Ht')}); 
G^(t,0 = -*([V'(i),V^t(i')]). (36) 

The first two Green's functions encode the dynamics 
whereas, the third, G^ (the Keldysh Green's function), 
encodes the evolution of the particle distribution. This, 
of course, is unnecessary in equilibrium, since the parti- 
cle distribution is known if the dynamic Green's functions 
are given. The two indices of the matrix G will be called 
"indices in Keldysh space" . 

Let us now discuss the specific propagators for the 
quantum dot Kondo model. One technical difficulty is 
that the Hamiltonian Hi contains, the spin operator S. 
To furnish a Wick's theorem, it proves convenient to 
factorize S in terms of canonical creation and annihi- 
lation operators,-. This is the "Abrikosov pseudo-fermion 
representation"E3 



(37) 



ct/3 



where f^a creates a pseudo-fermion with spin compo- 
nent a. To provide a faithful representation of the spin 
1/2 operator, we must impose the additional constraint 
Uf = 1. This is done—very conveniently by using the 
"Fedatov-Popov trick"c^E3, in which the f-electron field 
in the past is equilibrated with a heat bath at a com- 
plex chemical potential ^ = — i7rT/2. With this device, 
the contributions of doubly occupied and empty states 
exactly cancel. The occupancy of the f-electrons in the 
distant past is 



V- = fi-<7B + nrT/2) = 



1 



g-<TB/T+i7r/2 ^ I ■ 



where f{x) is the Fermi fimction. The magnetization of 
a free spin then recovers the Brillouin function 

M(°)(T,i3)= '^^/.-tanh(l). 



yielding a Curie susceptibility = dM/dB\B=o = y- 

The diagrammatic elements corresponding to the bare 
Green's functions are then 



m 

► — 

k,co 



(38) 



Gy(l,2) = -^(F^A(l)^t^(2)) =0, 



(34) 



for the conduction electron in the m-th lead and 
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G(u) 



CO 



(39) 



for the pseudofermion. Written out explicitly, 







(40) 



where eko- = Ck ^ crB and 





1 




27ri hfj5{ijj — \a) 
1 

CJ — Act —1(5 



(41) 



where we have introduced the notation Ao- — —aB. The 



functions hm{e) = h(t — Hm) and = h{\a 



-) give 



the occupancies of the band and f-electron states, where 
Me)^2/(e)-l = -tanhf-i;). (42) 



The quantity 5 is the usual infinitesimal regulator distin- 
guishing the retarded and advanced Green's functions. 



V. FEYNMAN RULES 

In developing the Feynman rules we need to consider 
two types of vertices. The "scattering vertex" generated 
by contractions of terms Hf = Hi{t+) — Hj(t-) in the 
expansion of the time ordered exponential is denoted by 






-<7ab ■ (^afj [{ — ) 



'■71 



LAA' 



(^)77'] • (43) 



(where the factor "z" derives from one factor of —i in 
the time ordered, exponential, and two factors of i trans- 
ferred from external legs). Here, a,/3,a,b refer to the 
spin indices and A, A', 7, 7' to the Keldysh indices of the 
outgoing and incoming propagators. The factors in bold 
are the matrices in Keldysh space that appear at the 
vertex. The shaded triangle denotes the vertex to which 
the Keldysh matrix ^ is to be applied. The Keldysh 
matrix entering into the vertex is actually symmetric un- 
der exchange of the Keldysh indices of the incoming, or 
outgoing fermions. 



(ti)aA' I77' 



+ (g) (ti)aa' = 

(ti)a7' «) I7A' + I7A' 



iri)xY (44) 



so that at each scattering vertex one has a choice as which 
pairs of incoming and outgoing fermion lines one applies 
the Keldysh matrices to. 
Diagrammatically, 



(45) 



The "measurement vertex" generated by contractions 
with Hf = ^[Hi{t+) + Hi{t^)] is denoted by 






<?ab ■ <?Q/3 [(^ 



lAA' 



( 2^)77' 



(46) 



(where the factor "z^ = —1" from the ti matrices is trans- 
ferred from the contractions of two external propagators.) 
Although the vertex (ti)aa' '9 (''"1)77' is not symmetric 
under interchange of incoming Keldysh indices, one can 
make it symmetric by adding in the additional "null" 
term Iaa' ® I77'. The expectation value of this vertex is 
zero, since it involves the product of two operators that 
are antisymmetric between the upper and lower Keldysh 
contour. By introducing this operator, using the identity 

(ti)aA'(8'(ti)^^'-|-1aA'(X'177' = (Ti)A7'<X'(-ri)T,A'-|-lA7'<X'l7A' 

and then removing the null operator on the right-hand 
side, one also has the freedom to interchange the Keldysh 
indices on the measurement vertex, i.e. 



(47) 



Our Feynman rules for the evaluation of the n-th order 
contribution to (Hi) are then 

• Construct a connected diagram with one measure- 
ment vertex ( ^ ) and n scattering vertices ( p3| ) 

• Associate factors g„j and G{uj) with the conduction 
and pseudofermion propagators. 

• Associate a factor — 1 with each conduction electron 
and pseudo- fermion loop. 

• Carry out the trace over spin and Keldysh indices. 

• Carry out the summation over internal momenta 
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• Carry out the integration over loop frequencies 
J Associate the convergence factor cosa;0+ = 
1 ^g»a;0+ ^g-»a;0+^ with propagators which form 
closed loops. 

• Associate a symmetry factor 1/n with diagrams 
with permutation symmetry group of dimension n. 

VI. SOME USEFUL DIAGRAMS 

A. Vanishing Diagrams 

In constructing the Feynman diagrams, there are a 
number of key diagrams which identically vanish. Any 
loop diagram which does not contain an insertion of the 
Keldysh matrix zti automatically vanishes, because it 
corresponds to the expectation value of a single, or a 
product of antisymmetric densities. For example, the 
simple fermion loop corresponds to the equal time expec- 
tation value of the difference between fermion number on 
the upper and lower Keldysh contour 



-i{nfa{t+) - nf„{t-)) = (48) 



In the frequency domain this result is obtained as follows 



/ 



d(jj r 

4n 



- — {—in + in) 
n 



0. 



(49) 



In this diagram, the first convergence factor catches the 
pole contribution from Ga, whereas the second conver- 
gence factor catches an equal yet opposite contribution 
from Gr. Two other examples of these cancellations are 



kma 




and 



= -i{nmka{t+) - rimkait-)) = 0. 



(50) 



0. 



where in this case, the cancellation does not require a 
convergence factor. 



B. Magnetization and Susceptibility 

The inclusion of the factor —iTi/2 into fermion loops 

gives the equilibrium expectation value of one particle 
quantities. For example, the equilibrium magnetization 
of the local moment is given by 



(7 

Y^''-h{K-inT/2) 



M„(T, B)=tanh( I 



Using the identity 
d 



(52) 



(53) 



we can differentiate the above result to compute the equi- 
librium finite-field susceptibility of the local moment 



_d_ 
dB 



3< 



dB 



^Tr f!^a,G{u;)a,G{u:] 



(54) 



In a similar fashion the magnetization of the electrons 

in the m-th lead is given by 

0-.=-'s/£-(ij.,g„,.,.,) 



2by: 



d.f„ 
dek 



= 2pB 



(55) 



This is the only case where we have to explicitly consider 
the effect of the magnetic field on the conduction electron 
lines. For all conduction electron propagators that do not 
loop back on themselves, we can carry out the summation 
over momentum k in advance of forming the Feynman 
diagram. We shall denote 



m 
► 

(0 



^ gm(fc, W) = ^m{'^)5o 



(56) 



(51) where 



gm(w) 



—inp 2niphm{oj) 
inp 



(57) 
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where we have used the large band-width approximation 

1 



E 



k 



(58) 



Note how the field dependence of the conduction electron 
propagator has been absorbed in the large band-width 
limit by a small shift in the conduction electron energies 
Cfeo-' — * Cfe- By making use of the identity 



In this approximation, the local magnetic susceptibility 
of the conduction electrons is given by 



1 /rfx^^M^^Va 
2./ x-e + iA 



£-iA 
2TTiT 



)-ln 



2-irT 



m,k(0 



(62) 



ak'co 



^Cr3gTO(w)(T3gn(wj 



5< >-/^ 

j di0[hm{u;) - hr,{w)] = 0. (59) ^fl(^) ^ ^ 



where A > 0, and tjjiz) = dlnr(z)/dz is the digamma 
function, we obtain 



P 
"2 



(a form of the Anderson-Clogston compensation theo- 
rem.) 



C. Kondo Polarization Bubbles 

Finally, for oiu' calculations, it proves useful to be able 
to compute the joint polarization bubbles between the 
f-electron and the conduction electrons. These bubbles 
control the renormalization of the interaction between 
the local moment and the leads. 

There are three basic bubbles: 



+ ■ 



2mT 



In 



27riT 



De~~' 
27rT 



2itT 



(63) 



The Keldysh loop is computed as follows 

dv 



27r 



-Tr 



— G^(W + I^)— gr„(l^) 



-iK „K 



-I 

= ^ ^ [1 - Kh^{tka')] 5{Xa -UJ- eka')- (64) 



dv 



SO that in the large band-width limit 



[1 - hcrhmiK - t^) 



(65) 



VII. DETAILED CALCULATION 



Note that the spin index a on each bubble refers to the 
spin component of the pseudo-fermion, and that each 
bubble is independent of the spin component {a') of the 
conduction electron propagator. The calculation of the 
first bubble proceeds as follows 



^Tr 
27r 



2,T\ 



_ f dv 
J Am 

-An 



'^a9m + ^o- .9: 
OJ + Cfe'cr' - Xa +i5 



(60) 



In the large large band-width approximation, this be- 
comes 



dv 



uj + v — \a +i5 



(61) 



We now discuss the evaluation of the interaction en- 
ergy, up to second-order perturbuation theory. The cor- 
responding Feynman diagrams are 



1. 



^ dJ' 




+- 



-{'Hl)ji,H = 



o:;:o +<::0::>+ 



4. 




= AF(i) + Af(2) + AF(3) + AF(4) 

where an additional factor of 1/2 in the last two terms 
comes from the coupling constant integration. 

In evaluating these diagrams, it is useful to divide the 
interaction vertex into an "Ising" and "x-y" component, 
writing 
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-Cafc • CTq/J — Jrr 



1 



(66) 

where cr^ = ^(cr^ ± icr^) are the raising and lowering 
operators. Notice that the amplitude for the spin flip 
terms is twice that for the Ising terms. Since the prop- 
agators are diagonal in the spin indices, diagrams (1-3) 
only involve the Ising component in the interaction. 

Taking account of the vanishing diagrams, we may 
compute the first three diagrams as follows: 



1. 




0(:> 

m _ 

{Jr + Jl)BMo 



9 2^ 4-^' 



2^4 

1 , - 



mni'Jnn 

{2pBY - 



dMo 
dB 



dB 



A = 



where the vanishing of the local conduction electron sus- 
ceptibility ( p9| ) causes the last term to vanish. 

The sum of the first three terms in the "free energy" 
is then 



(2) 



2^'^^^ -dB' 



(68) 




(70) 



These expressions can be divided into Ising (non-spin- 
fiip) and x-y (spin flip) components, AF(4-) = Ai^^^^^ + 

^F'^y where 



E Jmn f ^ I R K , K A \ 
A j rj \ ma n(7 ^ ma na j 





2 / J mn I 2y|. V n, — cj ' ma n, — cr J ^ 



(71) 

( The ^ prefactor in the first equation has been cancelled 
by the sum over the spin of the conduction electron lines, 
and the factor of (Jmn)^/4 is derived from the square of 
the Ising scattering amplitude. In the second equation, 
the spin of the conduction electrons is set by the spin of 
the f-electrons, so the ^ prefactor remains, however, the 
amplitude for two spin flips is ( Jm„)^, without a factor of 
1/4. ) Inserting the discrete expressions (^ ) and (|6|) 
we obtain 



mna ^ 

(67A^^± ^ Re ^ I 



(1 - Khnk){ha ~ 

Cfc - efe' 
(1 - h^„hnk){ha 



2aB 



(72) 



where Mo — tanh(-P) is the magnetization of the free 
impurity, and we have put Jr + Jl ~ J- We recognize 
these terms as the leading order expansion of the Free 
energy of a spin in a Weiss field -Be// = (1 ~ J)B: 

F* = -Tin [2cosh(B/r(l - J))] 

= -Tln[2cosh(i?/T)] + JBM^ - 77 (JS)'^ + ■ • ■ • 

z oB 

The Ising corrections to the energy are clearly indepen- 
dent of the voltage between the leads. 

The last diagram in the interaction energy can be ex- 
panded as follows: 

4. 



where we have denoted hmk = ^m(efc)- Following ear- 
lier discussion, the field dependence of the band-electron 
energies has been dropped in these expressions. 

As a simple check on these results, consider the zero- 
temperature, equilibrium limit of (72 ). In this limit, one 
can replace /i„i(efe') — sgnej;/, ~ a. The first term 
vanishes and the second term becomes 



AF± =Re V '^'»" [ (l-o-sgn(efc))(o- + sgn(efc')) 

ek-ek'-2<jB 

m.n.(7 ^ 



J' 



1 



£fc<0,efc/>0 



2B 



(73) 



(69). 

which is recognized as the second-order correction to the 
ground-state energy due to quantum spin fluctuations. 

To proceed further, it is convenient to introduce the 
notation 

Smna{x, B) — Re 27ri InF (x + 2<jB + /x,„ — /-t„) 
- (2; + 2aB + fim- Hn) In [ ) 



(74) 
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which has the property that 

^eiT^ai^ - Ain) = pS'„-^,^^{x, B) 

where the prime denotes derivative w.r.t. x. With this 
notation, using (63) and ( |65| ) we obtain 

^^(4) = E % / d^S'^nai^, B){1 - aMM^)) 

^^(4) =J2^ f d^S'^nA^, 0)(1 + aMM^)) (75) 

m,n 

where we have replaced Re/i_CT = —(tMq. Integrating this 
expression by parts, and using the fact that Smnai^D) = 
=FZ? we obtain 



2(P 
0(P) 



BMo In 



D 



2B 



(80) 



Again, the appearance of the equilibrium magnetization 
on the right-hand side of this expression has nothing to do 
with the value of the magnetization in the Hmit fj — > 0. 
With our method, we are unable to comment on the mag- 
netization in the limit J — *■ 0. However, the above ex- 
pression is very strongly suggestive that the limit J ^ 
returns to the equilibrium expression for the magnetiza- 
tion, in which case, 



M ^{1- J) tanh 
d 



(1 - J)B 



T 



^^(4) = - E 



■2 

m7i 



4J 



RL 



dB 



BMo In 



D 



2ttT 
BMo (in 





2B 


eV\ 






~Tj 


D 




/2B 


2ttT 


-4> 


\Y 



(76) +0{P). 



(81) 



There are two special limits of the above result. In the 
where the residual integral in AF^\^ vanishes under the limit T 0, we may take advantage of the asymptotic 



(4) 

spin summation. It is convenient to separate the logarith- form 
mic term away from the integral over Smnai'^)-, writing 

B 



In- 



where 



"^JmnBMo 



dx 



In 



D 
2^ 



2B fi^ 

T ' 



T 



D /2B V 

-lj2{2B + jV)\n 



(77) 

so that in this limit 



125 + 7^1 



27rDe 



(82) 



:.N9 77E27r«fTlnf(a; + crfe + 7u) (78) 
(2coshf )2 46 ^ / / V / 

^ fT,7 — ±1 

Combining the results (^, (||), and (|7|) our 

final result for the change in the Free energy AF — 
J '^{Hi) is then 



AF = -T(ln [2cosh(B(l - J)/T)] - F^qn - ^J^D 



M 



l-(Jfl + Jl)-2(J2 + J2)ln 



D 
2B 



In the zero field limit B 
$(?;) — 1 -|- 7, where 



0, 0(6, «) ^ 0(0, z;) 



(83) 



In 



^^JrlMoB 

'^{Jr + Jl)MoB 
0{P). 



D 



2ttT 
In 





2B 






~T 


~t) 


D 




f2B 


' 2tiT 







$(w) = Re 



dx 



^{v^x)-^(x) (84) 



(79) 



> 4cosh2(f) 

and i}{x) = V'(^ + *^)- ELud we have used the result 
dx 



Notice that the appearance of Feg„ on the right hand side 
of this expression does not rely on any assumption about 
the limiting value of the Free- energy or magnetization 
in the limit J ^ 0. The change in the magnetization is 
then given by 



4cosh^(f) 



^[x) = -(1+7) 



where 7 = 0.5772 . . . = —"0(1) is the Euler constant. In 
this limit, the magnetization takes the form 



M 



AM = (1 - J) tanh 



-4P — 



(1 - J)B 



B 
T 



T 



— tanh 



1 - 2{Jr -f Jl) - 4 {Jl + Jl) In 



2ttT 





/2B 




BMo 1 In ^ (t 






V 27rT ^ 







VlrI' 



In 



2ttT 



- $ 



(85) 
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VIII. DISCUSSION 

We should like to end with a discussion of some of the 
more controversial aspects of our work. There are two 
general topics that deserve discussion: 

* Does the magnetization of the quantum dot revert 
smoothly to the equilibrium value Mo = tanh (^) 
when the coupling to the leads is reduced to zero? 

* The issue of quantum coherence, and whether a 
current through the quantum dot dephases the 
Kondo effect. 

Both issues are areas where our current work differs in 
important ways with parallel work carried out by others 
in the field. 



1. Does the magnetization have a perturbative expansion? 

The first of the above items brings our discussion into 
contact up with the issue of whether a perturbative treat- 
ment of the magnetization in the quantum dot is valid 
out-of-equilibrium. Although strictly speaking, our ap- 
proach, based on the Onsager reciprocity relations, can 
not derive the leading term in the magnetization, the ap- 
pearance of a well-defined perturbation expansion, and 
the very form of our results 

AM = tanh(B(l - J)/T) - tanh(S/r) + 0(P) (86) 

suggests that the leading order term in the magnetization 
is just the Brillouin function tanh(i?/T). |— ■ 

Recent work by Parcollet and Hooley (PH)c3, Rosch, 
Passke, Kroha nd WolfleEd emd unpublished work by 
Kaminski and Glazman (KG)Eil casts doubt on the vi- 
ability of a perturbative approach. PH compute the lo- 
cal magnetization of a quantum dot by computing the 
steady-state Green function of the local moment, and 
they conclude that in the limit where the coupling of 
the local moment to the leads becomes infinitesimal, the 
magnetization does not revert to its thermal value, but 
instead acquires the limiting value 



^ {ip{b + v) + ^{b - v)) + 9^{b) 



(87) 



Mph{T,B,V) 
f{b,v) 

where 

and 9 = . A similar expression was derived by KG 

using a master equation approach. Results in the recent 
paper by Rosch et al. also reduce to the same form. 
This expression can not be expanded perturbatively in 



tanh I 



the coupling constants, and if correct, casts serious doubt 
in the whole viability of a perturbative approach to the 
non-equilibrium Kondo model. 

The KGPH result assumes that the relevant magneti- 
zation of the quantum dot is the polarization of the local 
dot spin, whereas we have argued that because the spin 
of the quantum dot partially delocalizes into the leads, 
one must consider the global magnetization. This leads 
to a fundamentally different philosophy about the way a 
quantum dot equilibrates its magnetization. Clearly, the 
local magnetization is not a conserved quantity, and it 
can be relaxed by its coupling to the leads. By exam- 
ining this problem, PH find a non-thermal polarization 
of the magnetic moment of the dot in the limit of zero 
coupling to the leads. We argue conversely, that if we fo- 
cus on the total magnetization then the coupling to the 
leads is not sufficient to relax to the appropriate steady 
state magnetization. In this case, a weak coupling be- 
tween the system and its surroundings is always required 
to relax the conserved total magnetization, which is oth- 
erwise conserved, even when the the dot and the leads 
are coupled. Let us call the coupling between the dot 
an its surroundings s. It is reasonable to suppose that 
physically reasonable results will be obtained so long as 
s is small compared with all physical energy scales in the 
quantum dot, such as the renormalized Kondo tempera- 
ture, but large compared with the adiabatic scale 6 which 
governs the rate at which interactions are turned on: 



Physical energy scales, (T^, eV, 



» s » 5 



In KGPH approach, in the limit of zero coupling to the 
leads, the local moment of the quantum dot is left in a 
non-thermal polarization. These results are most striking 
in the limit of zero temperature, when the polarization 
of the quantum dot is given by 



MpH = fJ-B X 



2-B(H 



eV/nB+2eB 
1 



{^iBB < eV/2) 
{iibB > eV/2) 



for fcsT << 2fiBB « eV, where we have restored 
V — > eV and B hbB. To preserve this non-thermal 
distribution, the spin must remain decoupled from any 
other thermal bath. In other words, the PH result cor- 
responds to the case where the coupling to an external 
thermal bath is zero, i.e s = 0. 

To get a feeling for the consequences of this result, 
consider a quantum dot at a bias voltage of lOOmF with 
9 — \. According to the PH result, the dot would have a 
magnetization M ^ O.l/is even in a field of 50 Tesla. In a 
real experiment, -this magnetization would relax back to 
one Bohr magneton, due to coupling to an external bath. 
In practice, such a coupling is always present, and can 
not be neglected when considering the expansion around 
the limit J = 0. Just as in equilibrium, providing these 
couplings are small compared with the relevant energy 
scales of the system- in this case the renormalized Kondo 
temperature- then a dynamical calculation of the mag- 
netization must include them. 
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It is difficult to reconcile the proposed failure of a 
perturbative expansion for the magnetization out-of- 
equilibrium with the existence of a perturbative expan- 
sion for the current and the total interaction energy of 
the dot. The current and interaction energy of the quan- 
tum dot appear to have a perfectly regular perturbation 
expansion. In the steady state, then the magnetization 
is determined from differential of the coupling constant 
integrated interaction energy: 



(88) 



While the Onsager reciprocity relation does not hold for 
charge current and interaction energy, we are able to re- 
late the current to the derivative of thefifflteraction energy 
with respect to the interaction energyO, 



dA 



(Hr 



(89) 



A=0 



where the above quantities are to be calculated by includ- 
ing a vector potential A into the scattering potential (but 
not into the measurement vertex of Hj), by replacing 



J, 



RL 



Jlr Jlrs 



Indeed, when we include these changes into our calcu- 
lation of Hj , we do recover the established second order 
expression for the current through the quantum dot (Ap- 
pendix D) 



3e^ 
h 



I=^{jRLfV. 



Thus a single energy functional appears to successfully 
yield both the magnetization and the current. From this 
point-of-view, a fundamental failure of perturbation the- 
ory would have to manifest itself in the interaction en- 
ergy. The smooth perturbative expansion of this quantity 
again suggests that the problems encountered by PH in 
their calculation are one of equilibration, rather than a 
break down of the perturbation theory at high tempera- 
tures or fields. 

^From the arguments of this section, it appears that 
the proposed failure of the perturbative expansion is tied 
up with two issues- the appropriate definition of the 
Kondo magnetization, and the the issue of how the mag- 
netization relaxes. If indeed, as we have argued, the total 
magnetization is the important variable for characteriz- 
ing the Kondo effect, then a coupling to an external bath 
is necessary in any dynamic treatment of the total mag- 
netization. It would be interesting in further work to 
examine whether such a coupling does dynamically re- 
lax the total magnetization to the energetically defined 
magnetization M. . 



2. Decoherence and Ward Identities 

Let us now turn to the second issue: that of decoher- 
ence. The general issue, whether a non-equilibrium quan- 
tum system loses its phase coherence is of great future 
importance, and will have bearing on non-equilibrium 
problems such as quantum computation. An issue here, 
fundamentally, is the question of whether driven non- 
equilibrium quantum systems can develop phase coherent 
structure, or not. 

In our quantum dot problem, we have used the On- 
sager reciprocity relations to argue that the energetically 
defined magnetization 



_d_ 



dJ' , 



is identically equal to the steady state magnetization of 
the quantum dot. 



M = M 



(90) 



We are able to clearly demonstrate a perturbative ex- 
pansion for Ai. Clearly even if the identity between the 
energetic and the physical magnetization were to fail, the 
perturbative expansion of A4 can still be used as a crite- 
rion for whether the Kondo effect enters a strong coupling 
regime at high voltages. The zero temperature, low field 
result 



M 



1- Jp- {Jpfln 



2B\{eV)\ 



(91) 



shows that the leading Kondo logs associated with in- 
tralead processes survive to arbitrarily high voltage. 

If decoherence is to prevent this divergence, at large 
voltages, then presumably, the above logs in A4 would 
have to be cut-off by the decoherence rate ~ I/e. 
This would require that inside the logarithms, we replace 

2B -> max(2B,?ir"i). 

For weak coupling to the leads, hr^^ ~ pJ^eV is the 
perturbative expansion for the current through the dot. 
Such cut-off effects will clearly appear if we extend the 
diagrammatics by dressing the fermion propagators using 
the Dyson equation 



(92) 



The main effect of the self-energy insertions into the 
fermion lines is to introduce a wavefunction renormal- 
ization of the magnetization spin vertex 



{Z ~ l)gp,B 



(93) 



However, in any controlled systematic treatment of the 
diagrammatics, it is important to also take into account 
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the corresponding renormalization of the magnetization 
vertex 



(94) 



Both of these diagrams are generated at the same time 
when one differentiates the interaction energy with re- 
spect to the external magnetic field. Now the important 
point is that these two vertices are related by a Ward 
identity associated with the conservation of total spin. 
In equilibrium, this leads to a Ward identity between the 
spin self-energy and the spin vertex 



ra;=o — 



di:{uj) 



doj 



leading to the precise cancellation of these two terms. 



+ <' 



duj 



= (95) 



a)=0 



This cancellation has important physical consequences. 
For example, this is the reason why the renormlized 
"meso-spin" S* (A) of the sort discussed in the introduc- 
tion 



S* 



preserves a zero relaxantion rate, whereas the bare spin 
S relaxes at a rate ^ J^ksT at high temperatures. 
The point is, that the magnetization of the local spin 
Miocai is rigorously given by 

Mlocal = TT-I - nj_ 

and this quantity therefore does not contain any vertex 
corrections. By contrast, the total magnetization enter- 
ing into the uniform susceptibility of a quantum dot does 
contain the above vertex. From this we learn that the 
rapid fluctuations of the local spin S described by the 
self-energy E are also responsible for the renormalization 
of the meso-spin S*{A) and its delocalization into the 
leads. This is why these two effects cancel. These con- 
siderations lead us to doubt any conclusions about the 
behavior of the D.C. biased quantum dot which rcsum 
the self-energy corrections to the pseudo-fermion lines 
without taking the corresponding vertex corrections into 
account. Any calculations based on the "non-crossing 
approximation" which assume that these vertices can be 
ignored will lose these subtle cancellations. 

It is useful to consider what will happen to the leading 
logarithm in the magnetic susceptibility if we continue 
to fourth order in the computation. The divergence we 
have found 



4(JI 



T 



In 



2ttT 



(96) 



corresponds to the diagram 




(97) 



where the logarithm derives from the f-conduction loops. 
When we go to fourth order, we dress the internal f- lines 
to produce the following three diagrams 



(98) 



The renormalization effects of this diagram will tend to 
cut-off the logarithmic divergence. However to the same 
order we also produce vertex counterterms 



AM, 



(4) 



(99) 



These vertex corrections correspond physically to the fact 
that the meso-spin S* (A) is now partially delocalized into 
the leads, with the rapid fluctuations of the original spin 
S now featuring as internal degrees of freedom inside the 
largely static meso-spin S*{A). It is our contention that 
the Ward identity associated with spin conservation will 
cause the decoherence effects in AM^^-* to be cancelled 
by the vertex renormalizations in AAfi , i.e. 



AM 



(4) 



amI^) = 0, 



and that furthermore, these cancellations will continue 
order by order. This is clearly an important point to 
confirm in future work. 

We would like to end by discussing the nature of the 
low temperature phase that develops in a symmetric 
quantum dot, where Jnfj = Jll- In our previous work, 
we suggested that the intersite coupling terms Jrl are 
marginally irrelevant at large voltage bias, so that the un- 
derlying low temperature physics of the large bias quan- 
tum dot is that of a two channel Kondo model. In such 
a picture, the magnetic susceptibility would have a loga- 
rithmic divergence of the form 



^oc— ln( ^ 



(100) 



where 



n/ieV). 



ji/ \cy J- More recently, Rosch et a]0 have 
suggested that the entry into the two channel Kondo 
regime will be interrupted at a temperature given by the 
decoherence scale 



Td 



V 



In' 
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In the former picture, a Curie-type magnetic susceptibil- 
ity would continue all the way down to temperatures of 
order T^, whereas in the latter picture, the magnetic sus- 
ceptibility would presumeably saturate at a much higher 
temperature. 

The survival of a Curie susceptibility at temperatures 
T «Td clearly relies on the existence of a renormalized 
composite "meso-spin" S* (A) and the vertex cancellation 
effects mentioned above. These same cancellation effects 
are intimately related to the ability of the meso-spin to 
transmit charge via internal fluctuations of the bare spin 
5, without the meso-spin undergoing fluctuations at all. 
A simple physical picture of this phenomenon involves 
equal but opposite currents of electrons and holes to flow- 
ing with the meso-spin. 




FIG. 7. Physical picture of how the "meso-spin" 5* (A) 
can carry charge from the left lead to the right lead, through 
internal motion of spin within its composite structure. When 
the dot-spin migrates into the leads, a charge current can 
flow without a spin current if the current of holes balances 
the current of electrons. 

At any given temperature, the original localized spin 
would then be in a state of rapid fluctuation, but the 
renormalized spin would remain essentially static until 
temperatures comparable with T^. 

3. Conclusion 

This paper has discussed the nature of the magnetiza- 
tion in the DC biased Kondo effect, highlighting it as an 
important central element in the debate about whether 
the DC biased quantum dot is decohered by a finite cur- 
rent, or whether it is able to enter a new highly non- 
equilibrium strong coupling state. We have proposed a 
slightly radical idea that the principle of virtual work can 
be used in the non-equilibrium steady state by consider- 
ing a large dot-lead system that is completely isolated, an 
idea that we showed is intimately related to our quantum 
reciprocity conjecture- the idea that subclass of variables 
in the non-equilibrium steady state continue to satisfy a 
reciprocity relation in time. Using these ideas, we com- 
puted an energetically defined magnetization and showed 
that it has a well-defined perturbation expansion in which 
the intra-lead Kondo logarithms are not cut-off by any fi- 
nite voltage. This result supports our earlier work, based 
on a point-splitting perturbation procedures. 

In the final discussion, we have attempted to bring out 
the most contraversial elements in the ongoing debate. 
We have argued that most of the proposed differences 
with our work stem ultimately from a failure to consider 
the total magnetization of the Kondo problem, rather 



than the local magnetization of the dot-spin. The total 
magnetization is not relaxed by coupling to the leads, and 
the physical susceptibility of real quantum dot always 
depends on a small coupling s to an external heat-bath. 
These difficulties can be avoided using an energetically 
defined definition of the magnetization, but they must be 
addressed directly if a dynamical calculation of the total 
magnetization is to be carried out. Finally, we have pro- 
posed that so-called de-coherence effects arise because of 
an incomplete and non-conserving treatment of the per- 
turbation theory. When spin conservation is taken into 
account, cancellations take place between the wavefunc- 
tion and vertex renormalization diagrams of the magne- 
tization operator that ultimately are associated with the 
delocalization of the dot-spin into the leads. We have 
proposed that these cancellations continue into the non- 
equilibrium state leading to a complete preservation of 
quantum coherence in the strongly biased quantum dot. 
These issues will clearly be the subject of future work. 

Experimentally, it is not yet feasible to measure the 
magnetization of individual quantum dots, however fu- 
ture improvements in the power of scanning Hall probes 
may make this a feasible experiment in the future. An- 
other approach to the experimental examination of these 
phenomenon may lie in a new configuration. Magnetic 
atoms inside a wire between two electron baths at dif- 
ferent voltages are exposed to two Fermi surface dicon- 
tinuities, and their low temperature physics is likely to 
be very ckisely associated with the D.C. biased quan- 
tum dot.LJ It may be possible by measuring the mag- 
netic susceptibility of pure silver wires containing rare- 
earth atoms with large Kondo temperatures, such as 
Cerium, to actually follow the magnetic aspects of the 
non-equilibrium Kondo problem. Experimental efforts in 
this direction would make an important counterparts to 
the theoretical debate we have outlined above. 
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IX. APPENDIX A: GREEN FUNCTIONS IN THE 
LARKIN OVCHINIKOV BASIS 

Wc begin by writing the Green function in the Larkin 
Ovchinikov Basis as 

G„^(l,2) = -i(Pvl/«(l)§^(2)) 

where 



*(1) = (Vs(1),^a(1)) 



(101) 



and 



- i(PVA(l)V^s(2)) = -- [(PV(l+)V'n2+)) 
+(PV(1-)V'^(2-)) + (PV(1-)V'^(2+)) + (PV(1+)V^(2-))] 

= [(Tv>(i)vn2)) - (Tv(i)vn2)) - imuHm] 

= Ga(1,2) 

With these manipulations, we see that we can write 



and 



m^^Hl)r.= [^^,f^). (102) 
Expanding this out gives 

*V(^V'a(1)V^a(2)) {PM^)i>^s{2))J ^'^''^ 

To transform this expression to its standard form, we 

need to expand the path-ordered Green- fimctions. Let us 
first carry this out for the (P'0s(l)V'^s(2)) term. Writing 



rn o^_ /^Gfl(l,2) Gk(1,2) 
^^^'^^"V Ga(1,2) 



X. APPENDIX B: FEYNMAN RULES FOR 
INTERACTION VERTICES 

To develop the Feynman rules for perturbation theory, 
we need to know the matrix elements of the interaction 
Hi the Larkin Ovchinkov basis. To compute these, let 
us consider the simplest case of an interaction 



^s(l) = ^(V(l+)+^(l-)) 

V'a(i) = ^ (V'(i+) - V'(i-)) 



Hi = ^j dld2V{l 



2)p(l)p(2). 



(105) 



(104) 



and then expanding the path-ordering along the Keldysh 
contout, we then have 



-i{PM^)^^sm = -- [(pv(i+)v^(2+)) 

+(P^(l_)Vt(2_)) + (PV>(l-)V^(2+)) + (PV(l+)V'^(2_ 

= -i [(Tv.(i)v^(2)) + ifm^H^)) + mi),i^H2)])] 

= -i{[m, V^(2)]) = Gk(1,2) 

in a similar fashion, we see the term (PV'a(1)V'^a(2)) 
term vanishes: 

-i{PM^)i^^A{2)) = -\ [(PV(i+)vn2+)) 

+ (PV'(l-)^t(2_)) - (PV'(1-)V^(2+)) - (PV(1+)V^(2_ 

= -\ [{Ti,{i)^\2)) + (Tv(i)vn2)) - {mi^\m] 

= 

By contrast, the diagonal terms are given by 

- i(PV5(l)V^A(2)) = [(P^(l + )^t(2^)) 

-(PV(1-)V^(2_)) + (PV(1-)V^(2+)) - (PV(1+)V^(2_))] 
= -|[(TV>(l)Vn2)) - (fV>(l)Vn2)) + ({^(l),^t(2)}) 

= -z(Ml),V^(2)Mti-t2) 
= Gk(1,2) 



where p{l) = '0^(1)V'(1) is the density. We can combine 
the density operator on the upper and lower contours of 
the Keldysh path to make a symmetric and antisymmet- 
ric combination 

P5(1) = ^(p(1+) + p(1-)) 

pa(i) = (p(i+)-p(1-)) 

i?he first combination is the measured value of the parti- 
cle number, whereas the second combination is the field 
that couples to an external chemical potential. In the 
Larkin Ovchinkov basis, these two operators can be writ- 
ten 

ps{l) = \ (V^s(I)V's(I) + V'^a(1)Va(1)) = i§(l)rif (1) 
PAiX) = (Vs(l)^^A(l) + ^t^(l)^s(l)) = *(1)*(1) 

))io that the "measurement vertex", ti/2 is off-diagonal 
whereas the response vertex is diagonal. This Keldysh 
structure holds for any one-particle operator. 

On the Keldysh contour we can define the measure- 
ment and response fields for Hi as follows 

ilf = \\ d\<mi(\ - 2)[p(l+)p(2+) +p(l_)p(2_)] 

Hf = lJ rfld21/(l-2)[p(l+)p(2+)-p(l_)p(2_)]. (106) 

The first of these operators is used to compute the expec- 
tation value of the interaction, but it is the the second 
operator which enters into the time-ordered exponential 
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and provides the scattering vertices that make up the XI. APPENDIX C: DERIVATION OF INTEGRAL 

Feynman diagrams. These can be re-written as IDENTITY 

Hf = 2 / 'ild2V{l - 2)[p^{l)p^{2) + -p"^(l)p'^(2)] In this section we prove two useful results. Firstly, that 



dld2V{l - 2)[p^(l)/(2) + p^(l)/(2)]. (107) 



m - 2 



2 2'KiT 



2 2'KiT 



(113) 
(114) 



The second term in involves the product of two an- 
tisymmetric operators, whose expectation value is zero. 
Since we only ever wish to calculate the expectation value 

of Hi, we can drop these terms. This enables us to write ^^^^^^ '/'(^) = d\nT{z)/dz is the digamma function and 
the interaction operators in the form second that 



1 



Hf = - dld2Viy^^^,{l - 2)i'A(l)*7(2)*7'(2)*A'(l) 



2 

i// = i / dld2y/^,^^^,(l-2)i'A(l)^7(2)*y(2)vI/v(l) (108) 



J 271 UJ - 

c 



c 



^( 



(115) 



where the Keldysh matrix elements for the measurement ^^ere $(cj) = D'^/{u}^ + D^) is a Lorentzian cut-off func- 



and scattering vertices are 



lAf,,w(l-2) = y(l-2) 



V3l,,^^,(l-2) = l-(l-2) 



'^1 ^ ■^i 
— (g) — 

2 2 



AA',77' 



„ 1 + 1 „ 

2 2 



(109) 



AA' ,77' 



In a Feynman diagram, each of these vertices is associ- 
ated with two propagators, introducing a factor (i^). The 
antisymmetric vertex also picks up a factor of — i derived 
from the path ordered exponential. The final results for 
the scattering vertices in the Feynman diagrams is then 



Vf,,^,^,(l-2) = y(l-2) 
VA>,^yil~2)^V{1^2) 



~2" ® T 



AA' ,77' 



-01 + 10- 



AA' ,77' 



These results are quickly generalized to spin dependent 
interactions. In this paper, we are interested in the 
Kondo interactions of the form 



tion. 



Consider the first integral /q. To evaluate it, we let 
^ and C both lie in the upper half complex plane. By 
completing the integral around the series of poles in the 
function f{z) at z = — «w„, we obtain 



lo = 2iT J2 



w„<0 

m 

n=l ' 



1 



1 



1 



1 



c 



(116) 



2 ' 2TTiT ' 2 ' 27riT / 

But the digamma function can be expanded as a series, 

1 1 



^iz)^-C + J2 



(110) 



r!=0 



n + 1 n + z 



(117) 



where —C — ip{l) is the Euler constant, so that 



°° / 1 



1 



n + y n + X 



(118) 



Hi = —Pa^apflSlpa^abi^b 

SO the interaction vertex is modified by replacing 

V{1 - 2) -> ^aab ■ Sap 

The corresponding interaction vertices are then 



(111) 



J 



— (g) — 

2 2 
iTi 
2 



1 + 1 



AA' ,77' 

in 

2 



(112) 



AA' ,77' 



where for clarity, we have omitted the indices from 
and V^. 



Substituting this into (116 ), we then obtain 



2 2TTiT 



2 2mT 



Let us turn to the second integral 



h 



did 



1 



(119) 



(120) 



where $((jj) = D"^ jiuP' + D"^). Expanding out the cut-off 
function, writing 



$(w) 



D 



1 



1 



uj — iD uj + iD 



then 
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h 



D f duj 
2i. ~ 



1 



1 



UJ — oj — iD 
in the limit where Z) 3> |Cli we have 



{D^-D} 



After making this gauge transform we obtain the follow- 
ing expression for the interaction energy to second order 
( at zero magnetic field B = 0), 



(- 



1 



1 



-) + {D^-D} 



But the second term is equal to the first one, because 

1 



(2) 



3J^ 



LR 



(k 

2^ 



/du) 
2^ 





1' 


1 


f duj 




1" 


~ 2 


Lu — iD J 


1 2^ 


^ 2 



uj + iD 



(the difference between these two integrals is the integral 
of an odd, with an even function, so it vanishes), so that 



TT 



/(-) - I 



1 



1 



Using Eq.(119), we get 



duj 
n 



m - 2 



uj ~ C, U! — iD 



(121) 



{[^f (e)7rf(e)+7rf(6)7rj^(6)] e"^ 

+ [^H(e)^f (^) e^} (126) 

where the factor 3 comes from summation of spins, 
Eq. (|66[) . By differentiating with respect to the vector 
potential, we obtain 

I''' = SJlnj / ^ {- [-f W K(^) ~ -He))] 

+ Kie){n^ie)-nt{e))]} (127) 



By using Eq.(|63|) and (pSj), the zero temperature current 
is then 



C ID 
^^^2 2'KiT^ ^ ^*'2 ^ 2^^ 



h J 2tt 



,1 



27riT 



,1 



27riT 



c 



27riT' 



In(^) 



Finally, since 



/ 



ZTT UJ — Q 



(122) 



(123) 



hence 



3e^^ 



(129) 



we reach 



/dhJ , , , , \ 1 
-/(.).>(.)— ^ = ^ 



^^2 27rir 



)-M;j^) 



(124) 



XII. APPENDIX D: DERIVATION OF CHARGE 
CURRENT AND CONDUCTANCE 

In this section, we derive an expression for charge cur- 
rent to the second order of J using the interaction energy 
(Hi). While the Onsager reciprocity relation does not 
hold for charge current and interaction energy, we are 
able to relate the current to the derivative of the inter- 
action energy with respect to the interaction energycSI, 



d_ 

dA 



(Hi). 



(125) 



A=0 



where the above quantities are to be calculated by includ- 
ing a vector potential A into the scattering potential (but 
not into the measurement vertex of Hj), by replacing 
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